Disorder and chain superconductivity in YBa 2 Cu307_,5 
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The effects of chain disorder on superconductivity in YBa2Cua07_i are discussed within the context 
of a proximity model. Chain disorder causes both pair-breaking and localization. The hybridization 
of chain and plane wavefunctions reduces the importance of localization, so that the transport 
anisotropy remains large in the presence of a finite fraction <5 of oxygen vacancies. Penetration depth 
and specific heat measurements probe the pair-breaking effects of chain disorder, and are discussed 
in detail at the level of the self-consistent T-matrix approximation. Quantitative agreement with 
these experiments is found when chain disorder is present. 
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Current understanding of the low energy electronic ex- 
citation spectrum of YBa2Cu307_5 (YBCO) is incom- 
plete. As with other high T c superconductors, YBCO is 
a layered compound in which the Cu02 layers are con- 
ducting, strongly correlated quasi-two-dimensional elec- 
tron gasses. The Cu02 layers have been extensively 
studiedEl in Bi2Sr2CaCu20s using angle resolved phcp 
toemission. While in-plane conductivity measurementsEI 
suggest strong similarities between the Cu02 planes in 
Bi 2 Sr 2 CaCu208 and YBCO, a complete description of 
YBCO is complicated by the presence of conducting one- 
dimensional CuO chain layers, about which relatively lit- 
tle is known. 

Since the CuO chains are unique to the YBaCuO fam- 
ily of superconductor, one expects to find their signature 
in a variety of experiments. Tunneling experiments,u 
for example, reveal a finite density of states (DOS) at 
the Fermi energy in the superconducting state, in con- 
trast to Bi2Sr2CaCu208 where a clear d- wave-like gap 
is observed.u Measurements of the penetration depth 
A c along the c-axis (perpendicular to the layers) find 
a power law dependence on temperature which is ma- 
terial dependent, Tu More dramatically, a well devel- 
oped pseudogapQ is found-, in the c-axis optical con- 
ductivity of YBa 2 C|U3|07_5El and YBa2Cu40sj3 but not 
in La2- :E Sr :E Cu04.H , E] In some cases, .such as in-plane 
anisotropic conductivity measurements £1 it is straightfor- 
ward to distinguish the contributions of the chains from 
those of the planes. In general, however, the influence 
of the chains is not trivial to understand. Calculations 
based on multi-band models suggest, for example, that 
interband transitions between the .plane and chain bands 
dominate the c-axis conductivity,!!!! and that the pseu- 
dogap seen in optical conductivity experiments reflects 
a shift in interband transition energies due to the open- 
ing of a gapli3 in the Cu02-layer Fermi surface. These 
predictions— are quite different from those of one-band 
models,E3~li3 in which disorder and inelastic scattering 
are assumed to play a key role in c-axis transport. For 
this reason, it is essential to develop a model which cor- 
rectly describes the low energy physics of the chain-plane 
system. 

Evidence concerning the chain electronic structure is 



indirect. In YBa2Cu307_,5, the chains are not continu- 
ous, but are broken into segments of finite length by the 
fraction S of vacant chain oxygen sites. In spite of this- 
there is a large anisotropy in the in-plane conduct ivitycl 
indicating that the chains are metallic. The absence 
of localization suggests that electronic states associated 
with the chains, while highly anisotropic, are not one- 
dimensionaL In the superconducting state, the in-plane 
anisotropytJ in the penetration depth A is nearly identi- 
cal in magnitude to the conductivity anisotropy, suggest- 
ing a significant supcrfluid density on the chain layer for 
temperatures T <C T c . What is truly remarkable, how- 
ever, is that the temperature dependence of the chain 
superfluid density — as measured in penetration depth 
experimental^ — is almost the same as that of the planes. 
The apparent similarity of the excitation spectra in the 
chain and plane layers is surprising given that the under- 
lying bands have completely different structures. 

Several different models have been proposed to de- 
scribe chain superconductivity. The simplest model con- 
sistent with d-wave chain superconductivity is the prox- 
imity model, in which the pairing interaction resides 
in the Cu02 planes and chain superconductivity occurs 
through the hybridization of plane and chain wavefunc- 
tions. The failure of this model (discussed below) to de- 
scribe penetration depth experiments,!^ has led several 
authorsEjOE!! to abandon the premise of a pairing in- 
teraction contained exclusively within the Cu02 planes. 

In this work, we show that a small amount of chain 
disorder is sufficient to reconcile the proximity model 
with experiments. We model a Cu02-CuO-Cu02 trilayer 
with a three-band tight-binding Hamiltonian in which 
the isolated chain and plane layers have one and two- 
dimensional dispersions respectively, coupled through 
single electron hopping. The resultant bands are three- 
dimensional hybrids of chain and plane states. Because 
the hybridization is weak for some values of the in-plane 
momentum k, localization effects are present, but are 
not sufficient to eliminate quasiparticle transport in the 
chains. We discuss disorder effects in chain superconduc- 
tivity in the context of two different experimental probes, 
of the low energy DOS: specific heat measurernents,Ej 
and penetration depth anisotropy measurements. 113 
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We consider a single trilayer with periodic boundary 
conditions along the c-axis. The annihilation operator for 
an electron in layer i with two-dimensional wavevector 
k and spin a is Ci k(T . The mean field Hamiltonian for 
such a model is (using the Nambu spinor notation) TL = 

E k c£H k c k with 
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amplitude t±_\ describes hopping between adjacent Cu02 
planes and tj_2 describes hopping between plane and 
chain layers. The dispersions £i k and ^2k describe the 
isolated plane and chain layers respectively, while A k is 
the superconducting order parameter for the plane layer, 
which is taken to have <i-wave symmetry. It has been 
suggested^ that a suitable tight binding model for the 
plane band is 

£ik = —2t± [cos k x a + cos k y a + 2t' cos k x a cos k y a 

+t" (cos2k x a + cos2k y a)] — fii (2) 

with t' = -0.2 and t" = 0.25, and a w 3l the lattice 
constant. The chain layer is modeled by 



^2k = —2*2 COS k y a — (12- 



(3) 



The parameters t\ and t 2 can be determined with 
some certainty from the magnitudes of the a (in- 
plane, perpendicular to chains) aad 6-axis (paral- 
lel to chains) penetration depths.tU The other pa- 
rameters are not as easily determined, but are con- 
strained by requiring that the Fermi surfaces be_ con- 
sistent with those of band-structure calculations S3 Ul- 
timately, our conclusions are not sensitive to the 
choice of parameters provided the above constraints are 
met. We take, as plausible: {t\, t 2 , Hi, ^2, t±i, £±2} = 
{60,400,-20,-600,40,60} meV. It should be empha- 
sized here that the chain bandwidth, 4*2, determined 
from Af,(T = 0) is veryjiear to that predicted by band- 
structure calculations ,lj consistent with the findings of 
positron annihilation studiesEil 

The order parameter is phenomcnological, with 



A, 



A (T)[g(k x ) - g(k y )}, 



(4) 



where the temperature dependence is given byc3 
A o (T)/A o (0) = tanh[T c A o (T)/TA o (0)]. Based on the 
location of the van Hove singularities in the tunneling 
DOS,clwe estimate Ao = 11 meV. Furthermore, we find 
that taking g(k) = cos(fca) — 0.3cos(3fca) gives approxi- 
mately the correct slopeO for A a (T) at low T. 

As we shall see, band structure plays a central role 
in determining the effects of chain disorder. The bands 




FIG. 1. Quasiparticle excitation energy along the nor- 
mal-state Fermi surfaces. The spectrum reflects both the 
d-wave structure of the order parameter, and the structure 
of the underlying band. 

(labeled a, b and c) have dispersions given by the positive 
eigenvalues of H k . In the normal state, e k = 6~> e k = e k ' 
and e£ = ejt where 
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and where £ k = £i k ± fj_i are the energies of bonding 
and antibonding combinations of the two planes. From 
Eq. ([|), it is clear that only the antibonding band mixes 
with the chain band, and that the bonding band (band 
a) is completely plane-like. In Fig. [I], we show the In- 
dependent gap in the superconducting excitation spec- 
trum, plotted along the Fermi surfaces of the normal 
state bands. Band a has the gap structure expected for a 
d-wave superconductor in a tetragonal system. Bands b 
and c have a more complicated spectrum, reflecting the 
structure of the underlying bands. 

It is simplest to start with a discussion of impurity 
effects in the normal state. Band a is unaffected by chain 
disorder, and we focus our attention on bands b and c. 
First, we emphasize that the behaviour of the two bands, 
given by Eq. (P), depends on the degree of chain-plane 
hybridization, as characterized by 



a k 



V2t 



±2 



16k - 67 



(6) 



In the limit a k 3> 1, the bands are degenerate relative 
to the coupling parameter t± 2 and e± « (^2k + £ k )/2 ± 
V / 2*_L2- This is the standard result for the level repulsion 
of a degenerate two-level system. The wavefunctions in 
this limit are even and odd combinations of the antibond- 
ing and chain wavefunctions, and the electron tunnels 
between the bands with a frequency 



a/r k = 2V2* X2 K > 1). 



(J) 



On time scales longer than r k , bands b and c are three- 
dimensional. In the limit a k <C 1, on the other hand, 
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e fc w min(£ 2 k,£ k ) _ °>k\/2t±2, e k « max(£ 2 k,£ k ) + 
«k\/2^±2 5 and the wavefunctions for bands b and c are 
predominantly plane or chain-like. The tail of band c 
near k x = 0, for example, is described by this limit. This 
weak mixing of chain and antibonding bands leads to a 
tunneling rate 

h/r k « V2t ±2 a k , (en, < 1) (8) 

which is much smaller than that of the degenerate case. 
In our trilayer model, is strongly k-dependent due 
to the different structures of the (one-dimensional) chain 
and (two-dimensional) plane dispersions. 

Because of the large variation in the tunneling rate 
ti/Tk, the effects of chain disorder are k-dependent. If 
the average length of unbroken chain is I, then we can 
estimate the scattering rate to be r 4 ~ = vf/1, where vf 
is the Fermi velocity in the chain layer. States for which 
the tunneling rate t^ 1 between layers is much larger than 
the impurity scattering rate r t ~ are considered three- 
dimensional, and are not susceptible to localization. On 
the other hand, chain- like states for which t^T 1 < r^ 1 are 
susceptible to localization. In our trilayer model, only a 
small fraction of chain electrons are localized by chain 
disorder. This observation provides a natural explana- 
tion for the large conductivity anisotropju seen in oxygen 
deficient YBCO. 

We now turn our attention to the superconducting 
state. In Fig. [l] we show the k-dependent supercon- 
ducting gap (defined as the energy required to excite 
a quasiparticle at the Fermi surface with wavevector k) 
produced by the pairing interaction in the Cu0 2 plane 
layers. The low energy tail in band c near k x = is of 
particular importance for the current discussion. States 
in this part of the band are predominantly chain- like, and 
the gap is 

£E»|A k |o£ K«l). (9) 

These low lying states have a significant influence on the 
low T properties of the system. In Fig. [| the clean- 
limit temperature dependent penetration depth has a 
pronounced upturn in AjJ~ 2 (T) at low temperatures. This 
reflects a sudden increase in the chain-layer superfluid 
density as T is lowered through ££. 

The r-kw temperature upturn in A^ is a generic 
featureEa-tj of proximity models of YBa 2 Cu307_,5, al- 
though the onset temperature for the upturn depends 
on the details of the model. In this work, we show 
that small amounts of chain disorder eliminate the up- 
turn. This is easily anticipated, since it is apparent that 
the states which have a small superconducting gap are 
chain-like, and are therefore strongly affected by chain 
disorder. The effects of disorder are two-fold: localiza- 
tion (discussed above) and pair-breaking. Pair-breaking 
occurs for chain-like states whose binding energy is 
less than the scattering rate h/T tr - Chain-like states are 
characterized by <C 1, and in this limit is al- 
ways smaller than the tunneling rate %/tu (assuming that 



A k < V^t±2). This means that any Cooper pair which 
satisfies the criterion r t ~ > t^ 1 for localization, is also 
broken by disorder. 

Localization effects will therefore be important in 
probes of quasiparticle transport. However, since the 
penetration depth is a measure of d.c. superfluid trans- 
port, the important effect of chain disorder is pair- 
breaking, which in this work is treated in the unitary 
limit of the self-consistent T-matrix approximation. The 
calculation is standard,Ej except for the fact that the im- 
purities reside only in the chain layer. In the Nambu 
notation, the self-energy S(w) is a sparse 6x6 ma- 
trix with nonzero elements connecting states in the chain 
layer only. The penetration depth then follows from 

A; 2 (0) - A" 2 (T) = £ 1 £ Tr[G(k, iu> n ) 

^ n k 

x 7/1 (k)G(k,iw n ) 7M (k)], (10) 

where Tr is the trace over the 6x6 matrix contained in 
the square brackets, G(k, iui n ) = [iuj n — H k — E(io; n )] , 
and 7/i(k) = h~ dUk/dk^. We emphasize that while 
this approach is reasonable for making predictions of su- 
perfluid and normal fluid densities, it cannot describe 
the enhanced-backscattering leading to localization of the 
normal fluid. EH 

In Fig. H we show the effect of disorder on A^ 2 (T) for 
different concentrations rii of unitary scatterers. We see 
that even small amounts of disorder are sufficient to elim- 
inate the upturn at low T. Chain disorder, on the other 
hand, has almost no effect on A~ 2 (T). In YBa 2 Cu307_a, 
the scattering results from a fraction 5 of chain oxygen 
sites which are vacant. Unfortunately, one cannot sim- 
ply equate 8 with the fractional impurity density rii intro- 
duced below because of the tendency of oxygen vacancies 
to cluster. Instead, it is more useful to consider the av- 
erage length of undamaged chain I (where I = a/rii) as 
the fundamental measure of chain disorder. 

A useful, and complementary, test of the model comes 
from examining the low T specific heat C v , which pro- 
vides a direct measure of the quasiparticle density. For 
d-wave superconductors, the electronic contribution to 
C v is well fitted by 

C v = l0 T + aT 2 (11) 

where 70 is proportional to the normal fluid density, 
and a depends on the k dependence of A k near the 
gap nodes. We evaluate 70 and a by fitting the DOS 
near the Fermi surface to N(u>) = Nq + Ni\u>\, and 
then evaluating 70 = NoN a k B Tr 2 /6 mJ/mol/K 2 , and 
a = 9YiY a fc|C(3) mJ/mol/K 3 , where N a is Avogadro's 
number, Ub is Boltzmann's constant, and £(n) is the Rie- 
mann ^-function. In the inset in Fig. ^ we show the effect 
of chain disorder on both 70 and a. We see that 70 in- 
creases rapidly with chain disorder, as a result of the 
increasing normal fluid density. These results are in ex- 
cellent quantitative agreement with experiments on high 



3 



U 0.5 




100.0 



T(K) 



FIG. 2. Temperature dependence of the penetration depth. 
Solid curves are X^ 2 (T) (upper curve) and \~ 2 (T) (lower 
curve) for ru = 0. Other curves are for A^~ 2 (T) for impu- 
rity concentrations m = 0.001, rii = 0.005, rii — 0.01, and 
rii = 0.02 from top to bottom. \~ 2 (T) is almost unchanged 
by chain disorder. Inset: Dependence of the coefficients of 
the low temperature specific heat on m. 



quality single crystals where 70 was found to be 1-3 
mJ/mol/K 2 for § between 0.01 and 0.05. 

Finally, we comment that we expect impurities to have 
little effect in c-axis transport, since the states which 
carry most of the current along the c-axis are those which 
are most three-dimensional, and therefore are least af- 
fected by chain disorder. In a series of conductivity 
experiments on Zn doped samples, Wang et alc3 have 
found that while the a-b anisotropy is strongly affected 
by small amounts of Zn, the c-axis conductivity is almost 
unchanged. 

In summary, we have shown that chain superconduc- 
tivity in YBa2Cu307_5 is well described by a proxim- 
ity model, if one accounts for the presence of disorder 
in the chain layer. Within such a model, hybridiza- 
tion of the chain and plane layers is strongly dependent 
on the in-plane momentum k, so that while most elec- 
tronic states are three-dimensional, a fraction of chain- 
like states are quasi-one-dimensional and susceptible to 
localization. We calculate the superfluid density in a self- 
consistent T-matrix approximation which captures the 
pair-breaking effects of disorder, aad find good agree^ 
ment with both penetration depthliZI and specific heatEHl 
measurements. 
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